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^ ■ Abstract 
J3 ■ 

We study generic solutions in a non-minimally coupled to gravity scalar field cosmology. It is 
shown that dynamics for both canonical and phantoms scalar fields with the potential can be re- 
duced to the dynamical system from which the exact forms for an equation of the state parameter is 
Q . derived. We have found the stationary solutions of the system and discussed their stability. Within 
^ ! the large class of admissible solutions we have found a non-degenerate critical points and we 
\ pointed out multiple attractor type of trajectory travelling in neighborhood of three critical points 
Ph| at which we have the radiation dominating universe, the barotropic matter dominating state and 
O ■ finally the deSitter attractor. We have demonstrated the stability of this trajectory which we call 
^ . the twister solution. Discovered evolutional path is only realized if there exist the non-minimal 
^ ! coupling constant. We have found simple duality relations between twister solutions in phantom 



and canonical scalar fields in the radiation domination phase. For the twister trajectory we have 
found an oscillating regime of approaching the deSitter attractor. 
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^ ! Observational evidences U 0] indicated that the current universe is dominated by an exotic 

Q\ \ form of matter with negative pressure called dark energy. The properties of dark energy component 
P ■ of the universe can be characterized by the equation of the state coefficient = p^/p^, where 
Pip and are pressure and energy density of dark energy component, respectively. In modern 
^ . cosmology the scalar field plays the major role in modelling of dark energy and in explanation 
of the accelerated expansion of the current universe (for review see isfl). In the simplest case the 
scalar field minimally coupled to gravity is assumed (so-called the quintessence idea [|J,|5i@,|T' 
In such a case the equation of the state parameter is > -1, but as it was noted by Caldwell 
the observational data also admit the possibility that < -1. Such a "phantom" form of a scalar 
field describing the dark energy component has many peculiar properties such as, for example, a 
big-rip singularity (energy density becomes infinite in finite time), the Lorentz invariance condition 



is then violated etc. [|9l.ll0l] 
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In this contribution we investigate dynamics of the cosmological model with the scalar field 
non-minimally coupled to the gravity with the positive and negative kinetic energy forms (i.e., 
canonical and phantom scalar fields) in the background of the flat Friedmann-Robertson-Walker 
(FRW) geometry. We point out interesting properties of a three phase model obtained within 
these class of solutions. They are interesting because they are generic and interpolate three physi- 
cally important phases of the evolution of the universe, namely, radiation, matter and dark energy 
domination in the evolution of the universe. Therefore, his solutions can be treated as a natural 
extension of the quartessence idea 111, 12, 13, 14]. 

In our investigations we apply the dynamical systems methods in exploring stationary states 



represented by critical points in the phase space as well as their stability [|15l . Il6l I17L 1 181 . Il9l |20| . 



21il22|,|23l|24|,|25|,|26D. We characterize all generic scenarios appearing in the case of the constant 
non-minimal coupling for both canonical and phantom scalar fields. In our dynamical study we 
relax the choice of the potential function. The presented approach to study the dynamics with 
the dynamical form of the equation of the state parameter is a different form the most popular 
one mainly used in the confrontation of the assumed model with dynamical dark energy with the 



observational data [|27l 12811 . While the authors who estimate parameters from the observational 
data postulate at the very beginning the form of the parameterization of the equation of state 
parameter w{z) as a function of the redshift z, in the presented approach such a form is directly 
derived from the closed dynamics of the FRW model filled by the non-minimally coupled scalar 
field. Moreover, basing on the twister solution one can derive approximated forms of the effective 
equation of the state parameter w(z) in three characteristic phases of the evolution of the universe, 
namely during the radiation, the baryonic matter and the dark energy domination. 

In the model under consideration we assume the spatially flat FRW universe filled with the 
non-minimally coupled scalar field and barotropic fluid with the equation of the state coefficient 
w„,. The action assumes following form 

s = ^Jd'x ^(^R - ^s'^^^M + ^R<f) - 2f/(0)j + (1) 

where = SnG, e = -1-1,-1 corresponds to canonical and phantom scalar field, respectively, the 
metric signature is (-, R = 6^^ -I- ^ j is the Ricci scalar, a is the scale factor and a dot 

denotes differentiation with respect to the cosmological time and U((p) is the scalar field potential 
function. Sm is the action for the barotropic matter part. 

The dynamical equation for the scalar field we can obtain from the variation 6S /6(p = 

+ 3H^ + ^R(p + sU'{(p) = 0, (2) 

and energy conservation condition from the variation 6S /6g = 

1 3 

fi = s-f + e3^H^(p^ + £3^H{4>^) + U{(p) + a» - -^H^- (3) 
2 

Then conservation conditions read 

3 9 

— W = p^+p,n, (4) 



H = 



^2 

y [(P^ + P^) + Pmi^ + Wm)] (5) 



where the energy density and the pressure of the scalar field are 
1 



1 



In what follows we introduce the energy phase space variables 



= y = — F — . z=—(p. 



which are suggested by the conservation condition 



-P0 + 



^0 + = 1 



or in terms of the newly introduced variables 

n^=y^ + e[(l - 6^)/ + e^x + zf] = l-Q„ 
The acceleration equation can be rewritten to the form 

'2 3 



H = - — [Pes + PeS) = (1 + Weff) 



2 

where the effective equation of the state parameter reads 
1 



Weff = 



1 - £0^(1 - D^jZ^L 

+(1 + w,„)(l - y") - el^(\ - 60z^ - 2^Ay\] 



where ^ = -^^^. 

K U((l>) d(p 



The dynamical system of the model under considerations is in the form [|16l] 

x' = -3x - n^z + s^Ay^[l - s6^z{x + z)] + e6^(l - 6^)xz^ + 

+^(x + 6^z)[e(l - 6^)(1 - w,„)x' + £2^(1 - 3w„,){x + zf + (1 + w,„)(l - /) 
y = -1^3;{^[l-£6^(l-6^)z2] + 6^/z}-el2^(l-6^)3;zH 



3 



+ '-y\ s{l - 6^)(1 - wjx^ + £2^(1 - 3w,n){x + zf + (1 + w,„)(l - /) 



z! = x[\ - £6^1 - 6^)z'], 

A' = -AHr-l)x[l-s6ai-6^)z% 



Table 1 : The location and eigenvalues of the critical points in twister quintessence scenario 





Weff 


location 


eigenvalues 


1. 


1 
3 




/i = -6^, /2 = 12^, /3 = 6^(l-3>vJ 


2. 


Wm 


x\ = 0,j* = 0,4 = 


/U = -1(1 - ± ^1 - f ^(H^), /2 = 1(1 + Wj 


3. 


-1 


x\ = ^Ay\f = 1,4 = 


/u = -i(3 ± V9 + e2a - 48^), h = -3(1 + w„,) 



where T = and prime denotes differentiation with respect to time t defined as 

A = ri_g6^(l_6^)z2l^. (17) 
dr J d In a 

If X is constant then we obtain the scaling potential exp (/l(/i) and the basic system reduces to the 
3-dimensional autonomous dynamical system in the case of the model with the barotropic matter. 
In the case without the matter the dynamical system is a 2-dimensional autonomous one. 
In the rest of the paper we will assume the following form of the function Y{X) 

= 1 - ;! (18) 

where a is an arbitrary constant and for which we can simply eliminate one of the variables namely 
z given by the relation 

where in the rest of the paper we take the integration constant as equal to zero. 

From equation (fTSi) and the definition of the function F we can simply calculate the form of 
the potential function 



f/(0) = f/oexp 



^l^^^+^e^j =f/oexp[-(<^ + r)2] (20) 



where yS is the integration constant. As we can see the dynamics of the model does not depend on 
the value of this parameter. In such a case we are exploring the solutions in the very rich family of 
potential functions. 

Following the Hartman-Grobman theorem ||29|] the system can be well approximated by the 
linear part of the system around a non-degenerate critical point. Then stability of the critical point 
is determined by eigenvalues of a linearization matrix only. In Table [H we have gathered criti- 
cal points appearing in twister scenario together with the eigenvalues of the linearization matrix 
calculated at those points. 

The critical point of a saddle type which represents the radiation dominated universe Weff = 1/3 
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is (xl = 0,yl = 0, (ip^ = ^) and the linearized solutions in the vicinity of this critical point are 

xiia) = ^-1— - i(l - 3w,„XAf - A\)]a-' + (1 - 3wj[xf + ^(Af - ^t)]ai-3"'-|21) 
yM = yfa\ (22) 
Ma) = ^i-2^-^|[xf-i(l-3wJ(4')-^*i)]a^ (23) 

In the case of canonical scalar field s = +1 this critical point exists only if ^ > and for the 
phantom scalar field e = -1 if ^ < 0. 

The matter dominated universe where Weff = i^m is represented by the critical point = 
0, y*2 = 0,A*2 = 0) which character depends on the value of the parameter d 

, ^ 16^ 1 - 3w„, 
d=l- 



3 (l->vJ' 

For <i > the critical point is of a saddle type and the linearized solutions are in the form 



.If 



r,-i(i-Wm) ( ^ ^ T ^ ^ T , r- 

X2{a) = |(1 + Vrf)[4'^ + -^(1 - w„){l - Vj)4'^]a-^(i-"''")^ - 



-(1 - yfd)[xf + ~{\ - wj(l + V^^^jfl^^i-"-")^!, (24) 



y2{a) = yfJ^'^'^-K (25) 

^ _ 2c.a-i(i--) |r ,) ^ 13 _ _ Vj)4)1«-|(i-.)V5 _ 
3(1-wJVjU' 2 J 

-[x® + ~(l - wj(l + Vj)4^]a^(i-^'"^^|. (26) 
For J < the critical point is of an unstable focus type 

- y/\d\xf cos (^(1 - w^) Vi^lnajj, (27) 

);2(a) = yffl^^^^^-^ (28) 

^2(a) = ^- ^-^a-^(^-"''"4[x« + - j(l - wJ4^1 sin(^(l - w^) Vl^lna) + 

3(l-wjVi^ 2 Q,4 2 J ^4 ; 

+^^(1 - wJ4^cos(^(1 -wj Vi^lnajj. (29) 

The final critical point represents the deSitter universe with Weff = - 1 is (Xj = 0, (jp^ = 1 , /I3 = 

0) its character depends on the value of A = 9 + s2a - 48^. 
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For A > the critical point is of a stable node type and the linearized solutions in the vicinity 
of this type critical points are 



x^ia) = ^^|(3 + Va)[x® + ^(3 - VA)Af]a-'^ - (3 - VA)[xf + ^(3 + VA)Af]a'^^) 
y,(a) = y; + (yf-y;)a-'^'--"'\ (31) 
Ua) = -^{[-f + ^(3 - VA)A«]a-^ - [xf + ^(3 + VA)i«]a^}. (32) 

For A < the critical point is of a stable focus type and the linearized solutions are 

[34'^ + i(9 + sa- 24^)Af] sin In a j + 4'^ Vi^l cos In a j|(33) 

(34) 



Xi(a) 


a 2 1 






= yl + (y 






Asia) 





|14 + 2^^3>in^— Inaj + cos ^— Inajj. (35) 

The solutions of the linearized system in the vicinity of each critical point Xi{a), yi{a) and 
Ai{a) can be used to constrain the model parameters thorough the cosmological data from various 
cosmological epochs. For example, the parameters for the solution describing radiation domi- 
nated universe (1) can be constrained from CMB data, and the solutions (3) describing the current 
accelerating expansion of the universe through the SN la data. Therefore one can estimate the 
parameters of the variability with redshift of true w(a) (see Fig©. It is possible because we have 
the linearization of exact formula in different epochs. 

The presented possibility of appearing twister type quintessence scenario is not restricted to 
the considered case of the Y{A) function (fTSl) . One can easily show that such a scenario will be 
always possible if only the following functions calculated at the critical points lllql 



fiA*) = (A*)\nA*) - 1) = const, = f'(A*) = const 

dA 

are finite. 

In this letter we pointed out the presence of the new interesting solution for the non-minimally 
coupled scalar field cosmology which we called the twister solution (because of the shape of the 
corresponding trajectory in the phase space, see Fig{T])- This type of the solution is very interesting 
because in the phase space it represents the 3-dimensional trajectory which interpolates different 
stages of evolution of the universe, namely, the radiation dominated, dust filled and accelerating 
universe. We found linearized solutions around all these intermediate phases as well as the solution 
of the separatrices which play an important role in the organization of the phase space structure. 
In is interesting that the presented structure of the phase space is allowed only for non-zero value 
of coupling constant, therefore it is a specific feature of the non-minimally coupled scalar field 
cosmology. 
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Figure 1 : Three-dimensional phase portrait of the dynamical system under consideration. Trajectories represent a 
twister type solution which interpolates between the radiation dominated universe (a saddle type critical point), the 
matter dominated universe (an unstable focus critical point) and the accelerating universe (a stable focus critical point). 




Figure 2: The evolution of Weff given by the relation (fT2l l for the non-minimally coupled canonical scalar field e = +1 
and the positive coupling constant ^. The sample trajectory used to plot this relation starts its evolution at tq - near 
the saddle type critical point (wefj = 1 /3) and then approaches an unstable focus critical point Weif - w,„ -0 and next 
escapes to the stable deSitter state with Weff = -1. The existence of a short time interval during which Weff - | is the 
effect of the nonzero coupling constant ^. 
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